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Q ■ Abstract 

CN ■ 

Wepointout that an erroneous derivation in the recent paper [Fetecauet al., Nonlinear Anal. RWA 12(2011) 1] yields 

^ . a correct solution by accident. Additionally, a number of misrepresentations and inaccuracies in the latter recent paper 

are identified, corrected and/or clarified in this Comment. Finally, a listing of recent papers in this journal that make 
Tij- . a mistake applying the Fourier sine transform, and thus present erroneous solutions, is given as an Appendix. 
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1. Introduction 



'^-H . The recent paper |ll|] considers the following partial differential equation: 

1-Q . which arises in the mathematical description of unidirectional planar motion of an Oldroyd-B (non-Newtonian) fluid 

'~7 ' resting upon an impermeable infinite plate situated at y = (see, e.g., \^ and the references therein for a complete 

derivation), with u - u(y, t) being the fluid's velocity component tangential to the plate. Here, v, A and A^ are constants. 

Furthermore, Eq. ([T]i is supplemented with the following initial and boundary conditions ol Eqs. (8)-(9)]: 
> ■ 

C^ ! du 

Tt . u{y, 0) = —{y, 0) = 0, y > 0; m(0, t) = t/(f), m ^ as y ^ oo, f > 0. (2) 

a\ ^ 

^^ , Here, we have omitted the condition du/dy ^ as y —> oo from yj, Eq. (9)] because it is extraneous ||2[]. 

t^^ ' The plate's velocity U(t) might appear to be an arbitrary function at this point, but it is not yfl. This is because, 

^^ , in following UJ, Eq. (8)], we have made a poor choice of notation in expressing the initial-boundary conditions. In 

fact, this has led to a significant number of errors when solving Eqs. ([T)-© by integral transform methods 12l|3l|4|]. 
Specifically, for Stokes' first problem, which is called the "Rayleigh-Stokes problem" in 111], the initial and boundary 
conditions are always "incompatible" at (y, t) - (0, 0) because of the requirement that the plate's motion be sudden 

k> \ (see the discussion in yfl). In other words, to ensure that the plate's velocity is identically zero prior to start-up but 

\^ • arbitrary (not identically zero) following start-up, as stipulated in Eq. (|2]i, one must write U(t) :- U(t)H(t), where 

. . ! U{t) is a smooth function that we are free to specify, and H(t) is the Heaviside unit step functionlJ This is not a "new" 

version of boundary condition, it is simply the mathematically-precise statement of the very same initial-boundary 
conditions for a suddenly moved plate stipulated in Eq. (|2]). Thus, for example, for Stokes' first problem, one takes 
U(t) = const., while for Stokes' second problem one may take U{t) - cos t. However, since H{t) is discontinuous at 
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f = 0, i.e., lim,^o+ H{t) + linif^o- H{t), it follows that U{t) is discontinuous at f = as well. Thus, great care must be 
taken in properly evaluating derivatives and integrals of the function C/(r) as written in Eq. ^. 

For what follows, it is important to realize that U{i) is a well-behaved function in physically-relevant problems. 
For example, we can take U{t) to be continuously differentiable for all f 6 R without loss of generality. Then, t/(0) 
is a well-defined quantity unlike U{0). Therefore, the claim that "for a greater generality we consider the boundary 
condition «(0, t) - U{t) with C/(0) = 0" 111 p. 3] is without merit when U{t) - const, because U{t) is discontinuous 
at f = 0. It appears that left-continuity of m(-, t) at the boundary is being implicitly enforced, without justification, in 
yj]. To the contrary, in the theory of non-smooth ordinary differential equations (notice the lack of smoothness on 
the right-hand side of Eq. (01 below), it is necessary to enforce r/^/if-continuity at f = (see, e.g., |^ p. 102] or ITl 
Chap. 15]), whence f/(0) ^ 0. 

2. How to correctly find the solution in the Fourier-Laplace transform domain 

Following [1], the next step in the analysis is applying the Fourier sine transform (in y with ^ as the "dummy" 
variable) to Eq. ([1]). Thus, it is claimed that the following ordinary differential equation (ODE) is obtained: 

^^+('^+c^f)^+yf-u.^^J-[yU(t) + aU'(t)], t>0, (3) 

or ot y 7T 

where u^i^, f) denotes the image of u{y, t) in the Fourier domain and a :— vA,-. Equation (O is subject to the trans- 
formed initial conditions u^i^, 0) = ~§f(^' 0) = 0. 

In order to be able to solve such an ODE, one must know a priori the degree of smoothness of the right-hand 
side. Does the classical solution of the ODE exist or should a solution be sought in the sense of distributions Jsl 
Chap. 2]? This is not made clear in Eq. (O above or in Ijj, Eq. (10)]. Instead, a cavalier application of the Laplace 
transform X (in t with q as the "dummy" variable and an over-line denoting the image of a function in the Laplace 
domain) is made in |ll]. Additionally, the authors implicitly use the relation X.{U'(t)}(q) - qU(q) - U(0) (incorrectly 
assuming, once again, that C/(0) - 0), which is only true for continuously differentiable functions |J8i §8.1, Property 
3]. Unfortunately, U{t) = U{t)H{t) fails to be continuous at r = 0, in general, as we explained above. Hence, the latter 
Laplace transform identity does not apply; moreover, the point value L'^(O) to be used in it is ill-defined for Stokes' 
first problem and related start-up problems. Consequently, the unjustified and erroneous assumption that U{Q) - 
renders the Fourier-Laplace domain solution of Eq. (O, i.e., |li Eq. (12)], incorrect. 

To remedy the situation, let us first note that the term U'{t), which comes form the mixed third-order derivative in 
Eq. ([TJ, constitutes a distributional derivative. Then, we have the following result: 

Lemma 1. U'{t) — U'{t)H{t) + U(0)6(t), where 6{t) is the Dirac delta distribution (generalized function). 

Proof. By the general result ISl §8.3, Eq. (4)] and the proper definition of the quantity U(t) from Sec. [1] U'{t) = 
[u(t)H(t)\ = U'(t)H(t) + lUMt), where [[[/]] := Um,^o+ U{t) - lim,^o- U(t). From the fact that U(t) is continuous 
at f = 0, we have lUJ = U(0) ■ 1 - C/(0) ■ = U{0). D 

By Lemma[ri upon applying the Fourier sine transform to Eq. ([T) subject to Stokes-type initial-boundary data, we 
obtain 

A^ + (1 + af)^ + vfu, =' ^ y ^ {vU(t)H(t) + a [u'(t)H(t) + U(0)6(t)]} , (4) 

rather than what is found in Eq. Q. Here, we have deliberately introduced the symbol = to stress that the equality is 

in the sense of distributions Jsl §2.3]. Note that U(0)6{t) = U(t)6{t) also ia §2.5, Example 2]; so, equivalently, we 
may write U'{t) - U'{t)H{t) + U(t)6(t), which is the result of formally differentiating U(t) = U{t)H{t). Now, to apply 
the Laplace transform to the right-hand side of the Eq. (|4]i, we first need to derive two identities. 



Lemma 2. £ [u'itWit)] (q) = qU{q) - U{<d). 



Proof. Following S §8.2], j^" U'{f)H{f)e-'i' At = ]^°° U'{t)&-'''At = J:[u'{t)]{q). Then, from S, §8.1, Property 3] 

and the fact that U{t) is a continuously differentiable function for all f e M (e.g., a polynomial or trigonometric 
function in Stokes-type problems), we have the desired result. D 

Lemma 3. Z {U'(f)} (q) = £ [u'(f)H{t) + U(Q)5{t)] (q) = qU{q). 

Proof. From |8i §8.1, Property 1], we have £ [u'{t)H{t) + C/(0)5(f)} iq) = £{u'it)H{t)] (q) + C/(0)X 16(1)} (q), which 
is just qU{q) - U{Q) + U{Q) ■ 1 = qUiq) thanks to Lemma|2]and S §8.2, Example 2(d)]. D 

Note that this result holds for any smooth U(t) and does not require the unjustifiable assumption that [/(O) = 
employed in yj] to "make things work out." In fact, we have shown the much stronger result that the value of U(t) at 
f = is irrelevant as long as U{Q) exists and is well-defined. 

Finally, upon acknowledging the initial conditions, we arrive at the correct solution of Eq. ([1]) in the Fourier- 
Laplace transform domain (with m7(^, q) denoting the image of u{y, t) there): 

K(^, q)^^J- , i^r^'^^Jl, ^ a ^(^)- (5) 

V n Aq'- + (1 + a^^)q + v^^ 

Though Eq. (|5]l is deceivingly similar to |jj, Eq. (12)], the former is the true Fourier-Laplace transform domain 
solution, while the latter is incorrect in general. Given the contradictory assumptions and lack of proof, it appears that 
111 Eq. (12)] is similar to Eq. (|5j only because the special case given in |4|, Eq. (3)] was "reverse engineered." 

3. Further deficiencies 

In yJ, Sec. 2], it is stated that "in order to solve a well-posed problem for such fluids one has to require an additional 
initial condition apart from the requirement that the fluid is initially at rest" Uj, p. 2] . This statement concerns the initial 
condition ^,{y, 0) = in Eq. (|2]l. However, for Stokes' first problem of the impulsively moved infinite plate, the fluid 
being at rest initially, i.e., u{y, t < 0") = (identically), trivially implies it has zero initial velocity u{y, 0) = 0, zero 
initial acceleration ^^(y, 0) = 0, zero initial jerk |jf (3', 0) = 0, zero initial jounce |jf Cv-O) = 0, and so on. Of course, 
these time derivatives are one-sided (with y > 0) because the solution itself possesses no time derivative across the 
plane f = as shown by the discussion above (see also JSl p. 577] and 121 p. 719]). None of these initial conditions are 
"additional assumptions" beyond the assumption of the fluid being initially at rest, which was made in the problem's 
formulation. Just because Eq. ([TJ is second-order in time and requires the use of two initial conditions stemming from 
the assumption that the fluid is initially at rest (rather than the single initial condition needed for, e.g., the second grade 
or viscous Newtonian fluids) does not mean that "additional assumptions" were made. 

Furthermore, an unreduced form of the yf-domain solution is presented in yJ, Eq. (17)], which appears to be able 
to give a complex-valued velocity if (1 -H a^'f' < 4-vA^^ in the expressions for ri,2. The concerned reader is referred to 
141, Eq. (7)], where it is shown that there are three cases to be considered and the solution is real-valued for all three. 

In yl Sec. 4.1], it is claimed by the authors that they "consider the case when the dimensionless relaxation and 
retardation times, A/t and A^/t - a/{vt) are much less than one" 111 p. 6]. However, this is a false statement. What is 
provided are results/or times that are long compared to both the relaxation and retardation time scales set by A and Ar, 
respectively. The dimensionless relaxation and retardation times, depending on the choice on non-dimensionalization 
scheme, could be given by, e.g.. A* - AUq/v and A* = Arlll/v, where U[) is a characteristic velocity. These are never 
introduced in yJ, Sec. 4.1]. 

In 111 Sec. 4.2], a scaling argument is given about the case of "/l/f, Ar/t » 1," leading to the claim that, in this 
limit, Eq. ([ll reduces to the analogous equation for a viscous Newtonian fluid if additionally A - Ar- However, it is 
clear from 141 Eq. (7)] that, when A = Ar, the solution to Stokes' first problem for the Oldroyd-B fluid is identical 
to the solution of Stokes' first problem for the viscous Newtonian fluid, independently of the sizes of A/t and Ar/t. 
The result in [;4] is mathematically stronger because it is a statement about the solution for A - Ar, regardless of any 
asymptotic approximation to Eq. ([T). The result is obvious once the solution to Eq. ^ is properly reduced, unlike the 
expression presented in |li Eq. (17)]. 



In addition, it is claimed that the correct solution for a second grade fluid executing the same motion (formally, the 
/I — > limit of the Oldroyd-B solution) "can be also obtained from [18, Eq. (2.5)] (integrating by parts the last term 
and taking V(t) - UH{t))" Ijj, p. 4]. Of course this is untrue as Hi Ref. 18] is known to contain erroneous derivations 
and solutions yD. Therefore, claiming ex post facto that a correct solution could be extracted from a wrong solution 
(once it is known what the correct solution is supposed to be) is fallacious reasoning. At any rate, the derivation of yl 
Ref. 18, Eq. (2.5)] suffers from the same logical inconsistencies explained in detail above. The reader is referred to 
^\M for more exposition on the issue. 

Finally, the list of references provided in jjj] fails to properly attribute previous results to the appropriate works in 
the literature. For the limiting case of a Maxwell fluid executing the same motion (formally, the /l^ — > limit in the 
Oldroyd-B solution), the energy analysis of Stokes' first problem is attributed to a 2007 paper yl Ref. 6] co-authored 
by two of the authors of 111]. The latter paper yl Ref. 6], the one under current discussion ll|] and a third paper by two 
of the authors of yj] in Int. J. Non-linear Mech. (vol. 44, pp. 862-864) in 2009, all neglect to mention that the energy 
analysis of Stokes' first problem for Maxwell fluids had already been completed in full in 2005 191]. 

4. Closure 

The failure to pose the initial-boundary conditions in Eq. (|2|l in a mathematically-precise form, combined with the 
inherently contradictory assumptions on the form and smoothness of the poorly-defined function U{t) from Eq. (|2]i, 
has lead to a large number of wrong solutions being published (see, e.g., ]l3i|4 i ll0ll for some proper corrections and 



also ]12I] for a longer listing of erroneous papers). Even a recent Erratum UllH claims that "for V{t) - V cos(ajO ... 
the solutions corresponding to the motion induced by an oscillating ... plate are recovered" Jlll p. 360] {V therein 
being U from the present Comment). Of course, limf^o+ cos(wO = 1 ^ 0, so the solution corresponding to the motion 
induced by a suddenly oscillated plate is not recovered. This is because of the lack of H{t) multiplying this expression, 
meaning its derivate calculated as in [|llL Eq. (5)] is generally incorrect (just as Eq. (O above is generally incorrect). 
Finally, it is also important to note that Eq. O is linear and has been solved (correctly) subject to the initial- 
boundary conditions in Eq. ^ as early as 1956 lll2ll (see also ]|2|, Sect. 2.2]). Therefore, it is unclear how the recent 
paper lUJ] u nder discussion here is "demonstrating the relevance and applicability of nonlinear techniques" (emphasis 
added)Til. 
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Appendix. Other erroneous recent papers on simple flows of non-Newtonian fluids in Nonlinear Anal. RWA 

Papers published in this journal that make the mistake of incorrectly differentiating the suddenly-moved plate's 
velocity (and therefore provide incorrect, unphysical solutions as discussed in i2ll3Ll4l ll0ll ') include, but are not limited 
to, the following: 

1 . F. Shen, W. Tan, Y. Zhao, T. Masuoka, The Rayleigh-Stokes problem for a heated generalized second grade 
fluid with fractional derivative model. Nonlinear Anal. RWA 7 (2006) 1072-1080. 

2. C. Fetecau, T. Hayat, C. Fetecau, N. Ali, Unsteady flow of a second grade fluid between two side walls perpen- 
dicular to a plate. Nonlinear Anal. RWA 9 (2008) 1236-1252. 

3. T. Hayat, C. Fetecau, M. Sajid, Analytic solution for MHD Transient rotating flow of a second grade fluid in a 
porous space, NonUnear Anal. RWA 9 (2008) 1619-1627. 

4. M. Khan, S. Wang, Flow of a generalized second-grade fluid between two side walls perpendicular to a plate 
with a fractional derivative model. Nonlinear Anal. RWA 10 (2008) 203-208. 

5. C. Xue, J. Nie, Exact solution of Stokes' first problem for heated generalized Burgers' fluid in a porous half- 
space. Nonlinear Anal. RWA 9 (2008) 1628-1637. 

6. M. Hussain, T. Hayat, S. Asghar, C. Fetecau, Oscillatory flows of second grade fluid in a porous space. Nonlin- 
ear Anal. RWA 10 (2008) 2403-2414. 

7. M. Khan, The Rayleigh-Stokes problem for an edge in a viscoelastic fluid with a fractional derivative model. 
Nonlinear Anal. RWA^IO (2008) 3190-3195. 

8. Y. Yao, Y Liu, Some unsteady flows of a second grade fluid over a plane wall. Nonlinear Anal. RWA 1 1 (2010) 
4442^450. 

The reader may verify that the above papers are erroneous by specializing the Fourier-Laplace domain representation 
of the governing equation given therein to the cases for which the correct form is known from i3[l4l ll0ll . then showing 
that the two disagree. Note that, at this time, only the paper in Item 5 above has been formally corrected thanks to the 
Comment paper of Jordan jlOll . 

What is more, the recent paper [L. Zheng, Y Liu, X. Zhang, Slip effects on MHD flow of a generalized Oldroyd- 
B fluid with fractional derivative, doi:10.1016/j.nonrwa.2011.02.016] purports to study once again a known linear 
problem. However, the article begins with the incorrect claim that "W = 0, pdV/dt - VT +pb" (instead of V • y = 0, 
pdV/dt = V ■ T +pb) are the "constitutive equations" (instead of the equations of conservation of mass and momentum) 
of an incompressible fluid. Skipping over the technical problems in between, consider the final conclusion of the 
paper: "the results indicated that the strongest shear stress occurs near the plate and the shear stress decreases rapidly 
with the increase of distance from the plate." Of course, this is nothing more than the imposed boundary conditions 
of a suddenly-moved plate and the decay of the velocity at infinity. It is unclear why the satisfaction of the imposed 
conditions can be the major conclusion of a scientific paper. 

So it goes. 



